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Abstract 

We prove an integrability criterion and a partial integrability criterion for 
homogeneous potentials of degree —1 which are invariant by rotation. We 
then apply it to the proof of the meromorphic non-integrability of the n body 
problem with Newtonian interaction in the plane on a surface of equation 
{H, C) = {Hq, Co) with {Hq, Cq) ^ (0, 0) where C is the angular momentum 
and H the energy, in the case where the n masses are equal. 
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1. Introduction 

Non-integrability of homogeneous potentials have already been a lot stud- 
ied mainly using Morales Ramis theory and its integrability conditions and 
Ziglin theory. These methods require a particular algebraic orbit of the cor- 
responding potential. With homogeneous potentials, there exist generically 
straight line orbits, corresponding to the Darboux points of the potentials. 

Definition 1. Let V : U C — > C be a meromorphic homogeneous 
function in gi, . . . , g„. We say that c G C" is a Darboux point if there exists 
a G C such that 

d 

-:—V(c) = aci Mi = 1 . . .n 
^q^ ^ ^ 
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We call a the multiplier, and we say that c is non degenerated if a 7^ 0. A 
Darboux point c is also call a central configuration in the case of the n body 
problem. 

These Darboux points correspond to homothetic orbits, which are explicit 
algebraic solutions of the differential equations 

d 

Qi = t:—V i = l...n 

With such orbits, it is already possible to prove some facts about non inte- 
grability, in particular in the case of homogeneous potentials. 

Theorem 1. (Morales, Ramis, Yoshida Q/,fl/,0/,Qy; // a meromorphic po- 
tential V is meromorphically integrable, then the neutral component of the 
Galois group of the variational equation near a particular (algebraic) orbit 
r is abelian at all order. Moreover, for V homogeneous of degree —1 and T 
a homothetic orbit associated to a Darboux point c with multiplier —1, the 
Galois group of the first order variational equation has an abelian neutral 
component if and only if 

Sp (VV(c)) C {^{k - l){k + 2), keN}. 

Here we want to study variational equations and their Galois group near 
another type of particular orbit that we often encounter when the potential 
is invariant by rotation. In particular, if there exists a plane of Darboux 
points, invariant by the rotation symmetry (this case in not rare), then we 
can build particular orbits with non zero angular momentum. Then we get a 
one parameter family of orbits for which we can apply Morales Ramis theory. 
For all of them, the identity component of the Galois group of the variational 
equation should be abelian, and then we can expect a much stronger inte- 
grability criterion than [ij. One difficulty is that the variational equation 
is too difficult to study in the general case, and then we will only make a 
complete analysis in the case we will call partially decoupled. We find very 
strong conditions, only two eigenvalues are possible instead of an infinity. 

The main results of this article are the following. 

We will first prove that if a homogeneous potential of degree —1 invariant 
by rotation is meromorphically integrable on a surface given by a fixed energy 
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and angular momentum, then the eigenvalue A of the Hessian matrix of a 
Darboux point with multiplier —1 should belong to the following table 



c 


A 


c = o 


A G {i(/c- l)(A; + 2), A; G N} 


C^E = -1/2 


A G {-P, keN} 


H = 


A G {i(A;- l)(A; + 2), keN\ 


{C, H) G C2 


AG {0,-1} 


(C,if) = (0,0) 


A G C 



where C is the fixed angular momentum and H the fixed energy. The com- 
plete statement is Theorem [6] with the table of Theorem [TJ in which there is 
an additional a priori hypothesis, the "decoupling condition" . Then we will 
apply this analysis to a well known case, the n body problem in the plane 
with equal masses. 

Theorem 2. The n body problem with equal masses is not meromorphically 
integrable on any hypersurface of the form C^H = a with a ^ fixed, nor 
on the hypersuface H = 0, nor on the hypersurface C = (H is the energy, 
and C the angular momentum). 

We see that this type of orbit allows us to study a new type of partial 
integrability: the case where the potential would be integrable only for a 
fixed value of the energy and angular momentum. This type of potential 
exists effectively as given in f lTOj) . This integrability table also gives indi- 
cations on which particular level we should focus (zero angular momentum, 
zero energy, and the case C^H = —1/2). This is for example helpful to do 
brutal search using Hietarinta [sj procedure for finding such potentials. The 
level C^H = —1/2 is also special because for example in the reduced 3 body 
problem we have an additional first integral in this level (the Jacobi inte- 
gral), although this first integral is not valid everywhere on C^H = —1/2. 
Theorem E] cannot solve all problems of this kind because of this "decou- 
pling condition" . A complete analysis of the 3 body case gives all the masses 
which satisfy this condition in Theorem [151 which are not always symmetric. 
A non integrability theorem like Theorem |2]is by the way immediate for these 
masses, except for (mi,m2,m3) = (1,5,1). 

Definition 2. We will call "norm" and scalar product the expressions 

n n 
||f ||^= ^^fj^ < V,W >= VjWi 

1=1 i=l 
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even for complex v,w (In particular, the "norm" can vanish for non zero v). 
We will say moreover that a matrix is orthonormal complex if its columns 
Xi, . . . , Xn are such that 

n n 

<Xi,xj>=j2mk{Xj)k = o vz,j \\Xi\\'=j2mi = '^ vi 
fc=i fc=i 

We note 0„ the complexified orthogonal group which is the group generated 
by these matrices, and SO„ the subgroup of 0„ of matrices with determinant 
1 (corresponding to rotations). In particular, the group 0„ conserve the 
"norm" . 

Definition 3. Let V be a homogeneous meromorphic potential of degree —1 
in dimension n>2. We note 

G^igeOn, V{g.x)^V{x)VxeC''} (1) 

We will call G the symmetry group of V. We will say that v e is in the 
equator of G if 

{ag.v, q; e C, g & G} 

contains at least a plane P of dimension 2 and that v e P. We will say 

that V is invariant by rotation if G contains at least a subgroup isomorphic 
to SO2. We will say that v is an eigenvector of G if for all g & G, v is an 
eigenvector of g. 

Theorem 3. Let V be a homogeneous potential of degree —1 in dimension 
n > 2 and G its symmetry group. Suppose there exists c in the equator 
of G such that c is a Darboux point of V with multiplier —1 and "norm" 
||c||^7^ 0. Then the variational equation near the conic orbit with parameters 
(C||c|p, £'||c||^) (angular momentum and energy) is given by 

t{-C^ + 2t + 2Ee)X + {-t + C^)X = i^-J) VV(c)i?0(t)X (2) 

where Reit) G SO2 with coefficients in C ^0, yj2E — C'^(j)~'^ + 20~^^ 

Proof. Let c be a Darboux point of V in the equator of G, with multiplier 
— 1 and "norm" 7 = ||c|p7^ 0. We note P a plane in G.c containing c. After 
rotation of the coordinates, we can suppose that c = (7, 0, . . . , 0) and that 
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P is generated by (7, 0, . . . , 0), (0, 7, 0, . . . , 0). A conic orbit for the Darboux 
point c corresponds to the orbit given by 

where ipt is given by 

Replacing this in the equation of energy conservation and angular momentum 
(the potential V restricted to the plane P is invariant by rotation), we get 

+ \i'4>'e' + ^ = - 

And after replacing, we get 

Knowing that the multiplier is —1, we have using Euler equation for V 

The variational equation is then of the form 

X = ^VV(i?,(,)c)X 

with R0(t) a rotation matrix. We also know that we are on some conic orbit 
(due to the fact that the homogeneity degree is —1) 

'^^'^ ^ 1 + e cos(^) 

with p and e some parameters depending on C, E. We get that cos(^), sin(^) 
are rational fractions in 0, <p. Then, with variable change — > t we get the 
following expression 

t(-C^ + 2t + 2Et^)X + (-t + C^)X = VV(i?0(t)c)X (4) 
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We know that the potential V is invariant by rotation. Then the matrix 
V^V^(-Re(t)c) equals to V^V^(c) after a basis change and gives 

V^ViReit)c) = R-^^V^V{c)Re(t) (5) 

Replacing this in (jlj) gives us the equation ([2]). □ 

Remark 1. The main difficulty of this variational equation is that it does 
not decouple after basis change. Indeed, we can make a basis change with 
some matrix P, but this matrix P should commute with the rotations Re{t)- 

Let us now give a proper definition of what will call integrable on some 
level of first integrals, as given in Theorem [2] 

Theorem 4. Let V he a homogeneous meromorphic potential of degree — 1 
in dimension n > 2. Let Ji, . . . , be meromorphic first integrals such that 

{J„J,} = V2,j 

where { , } is the Poisson bracket. We pose 

O = q) — > g{p, g, I{p, q)) with g holomorphic on 

the ring of holomorphic functions in p, q, J. We suppose that < Ji, . . . , > 
is a prime ideal on O and we pose K = Frac{0 / < Ji, . . . , 1^ >) the corre- 
sponding fraction field. Then the following functions are well defined 

• For all i = 1 . . . k , the functions ipi : K — )■ K, f — )■ {/, Jj}. 

• The function 



{n^i'i^)] f,g^{f,g} 



\i=l 



The functions ^ — )■ K which associate to fi, . . . , fn-k a sub deter- 
minant of size n X n of the Jacobian matrix (a matrix of size 2n x n) 
of h, ■ ■ ■ , h, fi, ■ ■ ■ , fn-k- 
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Proof. Let us write a representant of / G X as P/Q, P,Q E O. We just need 
to check that the value of the function does not depend on the choice of 
the representant. We consider hi, . . . , hk, gi, ■ ■ ■ , gk & O and we have 



Q + E 9sls 
s=l 



k 

-U k >| P+^hsIs ( k 

s=l J { s=l >! (Q + E 9sls)^ ^ 

s=l 

k 

k \-i P+Zhsh 

Q + 9sis {p, h) '-^ — {g, h) 

^ (g + E 9sisr 

s=l 



Q-\P,Ii} - PQ-'{Q,I,} = 



so the function is well defined on K. 

Let us consider /i,/2 € nf^^(y9~^(0) and we write P/Q a, representant of 
/i. Using the fact that {/j, = 0, we can do exactly the same calculations 
as before just replacing , Jj} by , /2}. Using the fact that the Poisson bracket 
is symmetric, we can do the same interverting the indices 1,2. So the function 
^ is well defined. 

Let us consider x one of the variables p,q, f & K and P/Q a representant. 
We have the classical formula 



So when we add to P and Q elements of < /i, . . . , /fc >, we are adding in the 
determinant a linear combination of . . . , dxlk- Using the fact that the 
determinant is multilinear, this will not change the value of the determinant 
because it contains the columns of the derivatives of /j, i = 1 . . . , /c. □ 

One need to be extremely cautious when manipulating these derivatives, 
because K is not a differential field, so we cannot conclude directly that all 
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notions we will need (Poisson brackets, independence) are well defined. For 
example, the Jacobian matrix itself we consider is not well defined on 
only its sub determinant of size n x n are. Remark that for the following, 
we will always consider a representant and forget the field K, and so this 
complicated definition will have no impact. This theorem is here to give a 
proper definition of integrability on some particular level of first integrals, 
and this complicated presentation has some advantages as it includes also all 
singular levels thanks to the prime ideal condition (we never ask for example 
that the first integrals be independent). 

Definition 4. Let be a homogeneous meromorphic potential of degree 

— 1 in dimension n > 2. Let Ji, . . . , 7^ be meromorphic first integrals such 
that {Ii,Ij} = Vi,j We say that V is meromorphically integrable on the 
manifold (/i, . . . , 1^) = if there exists functions Fi, . . . , Fn-k G K {K is 
defined as in Theorem H]) such that 

{H, F,} = oeK {/„ F,} = oeK {f„ f,} = o Vz,i 

and such that at least one of the sub determinants of size n x n of the 
Jacobian matrix of Ji, . . . , Ik, Fi, . . . , Fn-k is not in 7^ (this corresponds to 
the condition of independence almost everywhere). 

Theorem 5. Let V ^ be a homogeneous meromorphic potential of degree 

— 1 in dimension n >2. Suppose that an "angular momentum" 

C= ^ a.i,j{Piqj - VjQi) ttijEC 

k<i>j<0 

is a non trivial first integral ofV. Let us fix the value of energy H = Hq ^ 
and angular momentum C = Cq ^ 0. If V is integrable on this manifold of 
codimension 2, then V is integrable on the hypersurface C^H = CqHq. 

Proof. We consider the following transformation 

^ C'n _^ ^2n _^ (^^^ ^-2^) (g) 

We see that the transformation (f just multiply the Hamiltonian H — > a^H, 
and this does not change the integrability of H. Let us suppose that H be 
integrable on the manifold {H, C) = {Hq, Cq). We have 

Hivip, q)) = a'H C{v{p, q)) = a-'C {C'H){v{p, q)) = C'H 
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Then H is also integrable on the manifold {H, C) = {a'^Ho, a ^Cq). We also 
have 

U (a^Ho^a-^Co) = {{p,q) G C^" such that C{j,,qfH{p,q) = C^Ho} 

because CqHq ^ 0. This gives the theorem. □ 

Remark 2. We can see that the study of integrability on a specific manifold 
make sense only if this manifold is invariant by <y9, because if it is not the case, 
then our potential will be integrable on a manifold with higher dimension. 
Remark that the ideals < C — Co, H — Hq < C^H — CqHq > are always 
prime for CqHq ^ and V ^ 0, so integrability on these manifolds is well 
defined, contrary to the case < C^H > which will need to be splitted in two 
parts < C >,< H >. 



2. Integrability table 

Theorem 6. Let V be a homogeneous meromorphic potential of degree —1 
in dimension n > 2 and G its symmetry group. Suppose it exists c in the 
equator of G such that c is a Darboux point of V with multiplier —1 and 
||c|p7^ 0. We pose P the plane associated to c and G the group of rotations 
which stabilize the plane P. Let E be an eigenspace o/ V^V^(c) invariant by 
G. IfV is meromorphically integrable (respectively for some specific level of 
energy and angular momentum {H, G) ), then the following equation possess 
a Galois group whose identity component is abelian (respectively for some 
specific level of energy and angular momentum {H, G) ) 

t{-G^ + 2t + 2Ht^)X + {-t + G'^)X = XX H,G,XeC (7) 

where A is the eigenvalue ofV'^V{c) associated to the eigenspace E. 

Proof. This is a direct application of Theorem [31 We have a plane P and 
all vectors in this plane are Darboux points. The potential restricted to 
this plane is invariant by rotation (because the Darboux point c is in the 
equator). On the eigenspace E, the matrix V'^V{c) corresponds to Xld. We 
know moreover that the space E is invariant by the rotations Re(t) which 
corresponds to elements of G. We have then 



i?,-(J)VV(c)i?,(,) 



Xld 
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So the equation ([2]) on the eigenspace R simphfies and becomes equation 

m- □ 



Remark 3. The Theorem |6] has lots of hypotheses, but in fact only one of 
them is really restrictive. The existence of an invariant plane P such that 
the potential is invariant by rotation on this plane is common in practical 
cases. This often results by symmetry of the system. This is for example 
always the case in the n body problem. The restrictive condition is the 
existence of E invariant by rotation. In fact, this is a condition very similar 
to the codiagonalization constraint from Maciejewski, Przybylska they found 
studying potentials which are the sum of two homogeneous potentials. In 
fact, a potential invariant in dimension n by rotation can also be reduced to 
become a potential in dimension n — 1 which will be a sum of a homogeneous 
potential and the potential C'^/r'^. This new potential is not homogeneous 
and our condition correspond to the commutation of the Hessian matrices 
(at least on some non trivial subspace). 

Theorem 7. We consider the differential equation 

t{-C^ + 2t+{C^ -2)t^)X -{t-C^)X = XX C,XeC (8) 

This is a Fuchsian equation with 4 singularities, of Heun type f^J which 
corresponds to ([7]) with H = — 1. The corresponding Galois group is 
G = SL2{C) except for the following values of {G, X) 



G 


X 


G = 


Xe {l{k-l){k + 2), keN} 


G= 1 


X e {-k^, keN} 


G = V2 


X e {i(/c- l)(A; + 2), keN] 


G^{0,l,V2} 


A e {0,-1} 



Proof. We begin by the case G ^ {0, 1, \/2}. The equation (|8]) has exactly 4 
regular singularities on 

G' 

0,1,^^,00 
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We make Probenius expansion on these 4 singularities, and we find a loga- 
rithmic term ior t — and for i = oo. More precisely, we get 

2^. A- 2 - llA + 18 + 6C^ - 12C2 2 \ 

C2 (in t {^^e + Oit')^ - 2 - JjAt + o{e) 
^W-,(l + ^^ + 0(t-)) + 



These expansions are valid for A 7^ —1,0. In the case A = —1, we can 
compute explicitly the solutions and we find 

X{t) = Ci{t - C) + C2V(t-l)(tC2-2t + C2) 

The Galois group is then Z/2Z, abelian. In the case A = 0, we find the 
solution 

X{t) = ci + C2 (V(^-l)(*C''-2i + C2)) - 
/ (C--2). + l A 

The identity component of the Galois group is then C, abelian. Let us con- 
sider the case A 7^ —1, 0. Among the 3 solvable cases of Kovacic's algorithm, 
the only possible one with a logarithmic term is when there is a solution of 
the form 

X{t) ^^^p(^J F{s)ds^ F e C{t) 

If F has singularities of order more than 2 then X does not have a Puiseux 
expansion near this singularity. Impossible because all singularities are reg- 
ular. If the degree of F is positive, then the expansion at infinity is not a 
Puiseux series. Then the particular solution X{t) should be of the following 
form 

k 

x{t)^l[{t-ur' 
1=1 
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If rrii is not a non-negative integer, then tj is a singularity of X then ti equals 
to one of the singularities of the equation. This give even more constraints 
on the rrii because the Frobenuis exponents on 1, ^Sc^ 0, 1/2. On 0, the 
possible exponent is 2, and on infinity it is (the other ones correspond to 
the logarithmic behavior). This imply that the sum of the rrii be zero. But 
the rrii can never be negative, then they are all zero. The only left possibility 
is then X{t) = 1. We replace and we find A = 0, case already done. Then 
the Galois group is SL2{C). 

The cases C E {0, 1, \/2} correspond to confluences. These confluences 
are all regular (this has probably something to do with the fact that the 
system comes from a variational equation of a Hamiltonian system). The 
case C = has already been treated by Morales, Ramis, Yoshida [l|,[2j. Let 
us study the case C = \/2. This corresponds to the parabolic case (some 
study of this case has already been done by Tsygvintsev Q). The equation 
([8]) becomes 

2t{t - 1)X - (t - 2)X = \X 
There is a logarithmic term for the singularity t = 

x{t) = c,e [i + {^--\\^t + o{f)^ + 



for A 7^ 0, —1. In the cases A G {0, —1}, we find the solutions 

X{t) = ci + CsVt- 1(2 + t) X{t) = ci(t - 2) + CaVt - 1 

The Galois group is then Z/2Z in both cases, then abelian. We now look at 
the case A 7^ 0, —1. The possibles exponents are {2} at 0, {0, 1/2} at 1 and 



-- + -V8A + 9, -- - -V8A + 9 
4 4 4 4 

at 00. As before, we prove that we need a solution of the form 

k 



i=l 
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The possibles exponents outside infinity are always integers or half integers, 
all non-negative, then the sum of the rrii is a non-negative integer or half 
integer. Then 

4 4 2^ ^ 

We solve this equation and we find 

1 



A 



{k-l){k + 2) keW 



This is exactly the condition of Theorem [71 We now want to compute the 
Galois group for these remaining cases. We write the solutions of the equation 
in the following form (it is a hypergeometric equation, and the solutions can 
be written using hypergeometric series 2-^1) 



X{t) = ci 2F1 



C2 2^1 



1 1, 3 
-k, -k + - 

2 ' 2 2 



22 2 



, -t + 1 1 r+ 
t + i] Vt- It' 



These hypergeometric series are finite if the first bracket in 2-^1 contains a 
non-positive integer. For > 2, we see that either 1 — |fc or | — is a 
non-positive integer. Then one of the two functions is a polynomial. We 
always have a solution in C[t, y/t — 1], and then the identity component of 
the Galois group is abelian. 

Let us now study the case C = 1. The equation (j8]) becomes 

-t{t - lyx - (t - 1)X = XX 
The expansion on is the following 



X{t) = cit' 1 + A t + —A' — x + — ]r + o{t 
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6 / V 96 96 



C2 In t -A^ + -A + o{t^) - 2 - At 
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2 4 ' 



+ 



A 4 

and possess a logarithmic term for A 7^ 0,-1. The expansion at infinity is 

A , MA -5) , ^,^_3, 



Xit) = ci 1 



C2 In t A + 1 



A(A + 1) 
2t 



+ 0{t 



-2^ 



t + 1 



2t 12^2 
_ 4 + 11A + 3A2 
Tt 



0{t 

+ o{t-^ 
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Then there is always at least one logarithmic term for A 7^ —1. Remark 
that we already now that this equation has an abelian Galois group for 
A = 0,-1 (either using the limiting process of the generic solution for all 
angular momentum, or running Kovacic's algorithm for these specific cases). 
So now we will suppose that A 7^ 0,-1. We know that if the Galois group is 
not 5/72 (C), then it exists a solution of the form 



X{t) = exp F{s)ds^ F G 



The equation is Fuchsian and then X{t) can be written 



x{t) = i[{t-ur 



i=l 



The rui need to be non-negative integers except maybe at singularities. The 
exponents at 1 are +a/— A, —\/—\. Then one of the following equation is 
satisfied 



Then 



^+k = 



A 



or 



-A + A; = 



A; G N 



-(A; + 2)2 A;gN 

We add the cases A = 0, — 1 and this gives exactly the condition given by The- 
orem [71 We now need to compute the Galois group for these specific cases. 
We write the solutions of the equation in the following form (it is a hyper- 
geometric equation, and the solutions can be written using hypergeometric 
series) 



X{t) 



ci 2F1 

C2 2F1 



2 - iy/X, 1 - iVx] , [l - 2iVx 



1 



1 + i^A, 2 + i^/x] , \l + 2i\f\ 



t'(t-l)" 
1 -t) t^{t- 



+ 



These hypergeometric series are finite if the first bracket in the hypergeomet- 
ric series 2F1 contains a non-positive integer. We see that for A = — A;^ k G 
N*, it is the case for the solution in ci. There is always a polynomial solution 
and then the Galois group is always abelian. Remark that such a work can 
also be done using Kimura classification of hypergeometric equation which 
are solvable by quadrature in j^. 

□ 
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3. Algebraic potentials 

In the following sections, we will often need to consider algebraic po- 
tentials instead of meromorpliic ones. This is a problem because Theorem [1] 
deals only with meromorpliic potentials. This problem is often not addressed, 
except in Ziglin [ill, but in fact his procedure does not work. This is because 
making cuts in the complex plane does not allow after to make all possible 
monodromy paths. Then, the monodromy group will be reduced. It could 
have no consequences, but here there are important consequences because 
we absolutely need to be able to turn around the point in the variational 
equation (this is because for the two other singularities, the exponents are 
0, 1/2, and then if we restrict ourselves to these ones, the monodromy group 
will always be abelian). Let us now make a precise statement 

Definition 5. Let w be an algebraic "function". We define the critical set 

by 

= {x, w is not C°° on x} 

Theorem 8. Consider an algebraic potential V of degree —1 which is inte- 
grable with first integrals meromorphic in the positions, impulsions, and V, 
and a non degenerate Darhoux point c outside the critical set ofV. Then the 
identity component of the Galois group of the variational equation along the 
homothetic orbit associated to c is abelian at all order. It is also the case for 
an elliptic orbit if it exists under the conditions of Theorem O 

Proof. There are two important arguments to apply Morales Ramis theorem 

• The first integral need to have an expansion in series (or the quotient 
of two series in the meromorphic case) along the curve 

• The coefficients of this expansion will be functions of the time t, and 
the corresponding field will be the base field to be considered in Galois 
group computation. 

We have that V is homogeneous, then so is the critical set S(V^). If c ^ ^(y)^ 
then the whole orbit {c.(j),c.(j)) is not in ^{V), except maybe for = 0. 
It is also the case for an elliptic orbit using the notation ([3]) in proof of 
Theorem [31 So, a first integral meromorphic in the positions, impulsions, 
and V, is meromorphic everywhere on the particular orbit except possibly 
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for = 0. Let us prove now that the coefficients of an initial form of such 
ffist integral have coefficients meromorphic in 0, \/2E — C^0~^ + 20~^. 

Consider ffist the homothetic orbit O. To compute such an initial form, 
we need to derivate and evaluate the ffist integral / on O. We can write by 
definition 

-^(P) = /(P) Q') ^(Q')) with / meromorphic 

The problem is of course the derivation in V on O. We know it is deriv- 
able everywhere on O except maybe for = 0. But we also have that the 
derivatives of V are homogeneous functions. So we have 

5V(c.0(t)) = Hif)-^-^ (3eC 

where correspond to a derivation of order k (and /3 depend of course of 
the chosen derivation). We see that even if there is a singularity on = 0, 
the function is still meromorphic on = 0. The coefficient (3 is well de- 
fined because the derivative is well defined for 7^ 0. In the case of elliptic 
orbits, the potential is invariant by rotation on the plane generated by the 
elliptic orbit. So the only problematic point is also (0, 0). As given in Theo- 
rem [3|, the functions cos(^), sin(^) of the angle are functions meromorphic in 
0, ^2E - C^cj)-^ + 20-1 and then, so are the coefficients of the initial form. 

To conclude, we now need to precise that, as Morales and Ramis using 
Kimura table |8| have done, we are in fact computing Galois group over 
the base field C(0), not on the field of meromorphic functions. This is not 
a problem here because in all cases, the variational equation is regular at 
infinity, and then so are the first integral. We know that a meromorphic 
function on C is in fact rational, so the coefficients of the initial form are in 
fact in 

C (0, ^/2E^^C^^fF^^2^'^ 

This is just an extension of degree 2 of the field C(0), so the identity com- 
ponent of the Galois group will be the same. □ 

4. The case of dimension 3 

In the particular case of dimension 3, we get 

Theorem 9. Let V{x,y,z) be a potential meromorphic in a/x^ + y^, 

a/x^ + + 2;^, 2^ and homogeneous of degree —1 in dimension 3 (this imply 
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that the symmetry group ofV contains Z/2Zx 02^. Suppose that V{1, 0, 0) ^ 
0,00. If V is meromorphically integrable, then it belongs to one of the fol- 
lowing families 

V = " aeC* V = ^ b G C* (9) 

This theorem is almost the best we can have (with a reasonable state- 
ment). To apply our previous theory, we need an invariant plane on which 
the potential is invariant by rotation. Such invariant plane comes here from 
the symmetry in z. The constraint V^(1,0,0) 7^ 00 cannot be removed, but 
the constraint V{1, 0, 0) 7^ could maybe be removed with a lot of additional 
work. There are two keys which allow us to give such a complete statement, 
which are the fact that the decoupling condition is always satisfied, and then 
that the potential can be reduced on a plane for which an almost complete 



classification is already done in 10| (for a finite number of eigenvalues). 



Proof. The potential V possess a symmetry group G such that 
G D 



cos{9) 


—sin{9) 









sin{9) 


cos{e) 


Dil 


1 


















We consider P the plane z = 0. This is an invariant plane because dzV{x, y, 0) 
= thanks to parity in z. Using the hypotheses, the restriction of V to 
the plane P is not zero or infinite. The point c = (1,0,0) is then a non 
degenerated Darboux point, and c is in the equator of G because c E P 
and G.c D P. The matrix V^V^(c) contains a stable subspace of dimension 
2 associated to P. Then the supplementary space generated by the vector 
(0, 0, 1) is also an eigenspace. The rotation group generated by the rotations 
around the 2;-axis let invariant the vector (0, 0, 1). The conditions of Theorem 
|6]are satisfied and the "vertical" (normal to the plane P) variational equation 
is then 

t(-C2 + 2t + 2Et^)X + (-t + C^)X = d,,V{c)X 
This equation is integrable for all angular momentum G only if 

9,,V(c)G {0,-1} 

with c with multiplier equal to —1. We now restrict our potential to the 
plane P : y = 0. The potential V is invariant by rotation around the 2;-axis, 
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then P is an invariant plane and we consider the restriction V : P ^ <C. 
The restriction of the function -y/xMT^ to y = gives a bivaluated function 
whose values are +x or —x. We can choose one or another for the restriction 
V (because V should be integrable for both possibilities anyway), and we 
choose arbitrary 

y/x^ + \y=o= X 

The function V{x, z) is then meromorphic in x, y/x"^ + z^, z. It possess a 



Darboux point c = (1,0), and it is non degenerated. Using [lOj, and then 
applying the symmetry group, we find that if V is meromorphically inte- 
grable, then it should be of the form ([9]). These potentials effectively possess 
an additional first integral, / = p^x — PxZ and I = Pz respectively, and they 
are independent with energy and angular momentum almost everywhere. □ 

Remark 4. We can see here the importance of the symmetry group struc- 
ture in the study of integrability. Here, the vertical variational equation is 
simple because in dimension 3, a group of rotations (except SO3) always 
possess an eigenvector. This is no more the case in dimension 4 and higher. 
In particular, the complexity of the variational equation is very linked to the 
symmetry group, and if it is too complicated, we will need additional prop- 
erties on the matrix V^\^(c). As we will see after, in the n body problem, an 
explicit decoupling condition appear because the symmetry group contains 
the rotations 

Jl - [ cos6'/„ -sin6'/„ 
^ \^ sin6'/„ cos 6 In 

and this group do not possess a common eigenvector of eigenvalue 1. 

Now we can ask if the Theorems [6H71 are really "useful" and not only 
purely theoretical possibilities with no examples. Does exist effectively some 
potentials that would be integrable only for a specific value of energy and 
angular momentum? Using Hietarinta |5| direct method and then our non 
integrability approach, we find the following potentials 

Theorem 10. We consider the potentials 

^ + ^ x^ + y^ + z^ 

The potential V\ is integrable for zero angular momentum C = 0, but not on 
any other hypersurface C^H = a, a G C* (the question about integrability 
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on H = is still open). 

The potential V2 is integrable on the hypersurface H = of zero energy, but 
not on any other hypersurface C^H = a, a G C*, nor on the hypersurface 
C = of zero angular momentum. 

Proof. The first integral of Vi is given by 

(xp^ + ypy)pz z 



The potential Vi possess a Darboux point (1, 0, 0), and the associated eigen- 
value is A = 2. Using integrability table of Theorem[71 we have that this value 
is only possible for the hypersurfaces C = and H = 0. Then Vi is not mero- 
morphically integrable on any hypersurface of the form C^H = a, a G C*. 
The first integral of V2 is given by 

h = {x^ + - z^YpI - 4:z{x^ + - z^)pz{xpr, + ypy) + Az^{xp:, + ypyf 

The potential V2 possess a Darboux point (1, 0, 0), and the associated eigen- 
value is A = 2. Using integrability table of Theorem [71 we have that this 
value is only possible for the hypersurfaces C = and H = 0. We know it 
is integrable for H = 0. Suppose it is integrable for C = 0. Then we could 
reduce the potential by rotation and we would obtain the following potential 

f. ?! + 
V2 — 3 

This potential possess a Darboux point (1, 0) and the associated eigenvalue is 
A = 2. But in this case, it already has been proved in [lO|] that the potential 
should belong to one of the following families (after rotation) 

V=^+^ a,beC* ^ = ^M±4+^^ aGC*,e = ±l (11) 

qi q2 [qi + qi + m2 

The second case is impossible because it is always complex. For the first one, 
we apply a rotation to V2 of angle 6 

y ql + ql 

(cos(6')gi + sin(6')g2)^ 

and this never coincide with expression (ITTi) . Then V2 is not integrable on 
the hypersurface C = 0. □ 
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5. Application to the n body problem 

We consider V the potential of the n body problem in the plane 

^^ym^ (12) 

with positive masses mj, qi G C^. The symmetry group is (at least) 

cos 9 In — sin 9In 



sin 6 In cos 6 In ' ' ^ ^ 

Let c be a Darboux point with multiplier —1 and such that ||cp7^ 0. Then 
c is in the equator of G and we can build a conic orbit (by definition, the 
mutual distances between the bodies are not zero). For the following, we will 
pose 

using notation mj+„ = mj. Remark for the following that the potential of 
the n body problem as given by (fT2|) is not reduced at all. This means in 
particular that the kinetic part is 



2m,- 



and so does not correspond exactly to the case we studied before. Still, it 
is almost the same and we just have to make a variable change like Pi — > 
Vi\p^i- The matrix V^l^(c) becomes in particular the matrix given by (fT3ll . 



5.1. General properties 

Definition 6. Let V be the potential of the n body problem with positive 
masses mj, c a Darboux point with multiplier —1. We will say that the 
variational equation near a conic orbit is partially decoupled if it exists a 
non trivial vector space V and A G C such that 

Wv = \v Vt; G ■(> 

and V is stable by the rotations 

^ _ / cos Qln — sin Qlr, 
^ sin Qlr, cos Qlr, 
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Remark 5. This definition exactly corresponds to the existence of a non 
trivial eigenspace satisfying Theorem |H1 



Theorem 11. Let V be the potential of the n body problem with positive 
masses uii, c a Darboux point with multiplier —1 and W G M2n(C) the 
associated matrix (given by equation (|T3|l ). The variational equation near a 
conic orbit is partially decoupled if and only if it exists a vector v G C^" \ {0} 
and A G C such that 

Wv = J-^WJv = Xv (14) 
where J G M2n(C) is matrix of the canonical symplectic form. 

Proof. Suppose at first that v is not an eigenvector of Rg (these matrices com- 
mute so they have the same eigenvectors). We just have to take V = {v, Jv) 
because the space generated by Rgv, \/9 is a 2-dimensional space which con- 
tains {v, Jv) {V is always 2-dimensional because v is not an eigenvector of 
J = Rt,/2)- Using the hypotheses, v and Jv are eigenvectors of W with the 
same eigenvalue, so V is an eigenspace of W stable by the rotations Rq. liv 
is an eigenvector of Re, then we take V = C.v and V is an eigenspace of W 
stable by the rotations Rg. 

Conversely, if we have an eigenspace V stable by the rotations Rg, we take 
any vector v E V and it satisfy (IT^ because J = Rn/2 and then Jv G V, and 
so it is also an eigenvector of eigenvalue A. □ 

Theorem 12. Let V be the potential of the n body problem with positive 
masses rrii, c a Darboux point with multiplier —1 and W G M2„(C) the 
associated matrix (given by equation (fT3|) ). If the variational equation near 
a conic orbit is partially decoupled then the matrix W of (1131) has a double 
eigenvalue. 

Proof. If dim{V) > 2 then by definition the matrix W has a double eigen- 
value. Let us consider the case dim{V) = 1. The corresponding vector have 
to be a common eigenvector of J and W. The eigenvectors J are of the 
form {w,iw), w G C". In particular, they have zero "norm". But if W has 
only simple eigenvalues, then W is diagonalizable and using Theorem 6 of 
Til , W is then diagonalizable in an "orthonormal" basis. So if W has an 



eigenvector with zero "norm", then this eigenvector is a linear combination 
of two eigenvectors and this implies an eigenspace of dimension greater than 
2 and then a double eigenvalue. □ 
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Theorem 13. Let V be the potential of the n body problem with positive 
masses mi, c a Darboux point with multiplier —1 such that the bodies are 
aligned and W G M2n(C) the associated matrix (given by equation f lT3|) ). 
We suppose that W is diagonalizable. Then the variational equation near a 
conic orbit has a Galois group G such that 

G^G with GcCx 5p(2)"-2 
where Sp{2) is the 4 dimensional symplectic group. 

Proof. For an aligned Darboux point, we have the following property (found 
by direct computation) 

Then W and J~^WJ commute. Then it exists a common eigenvector basis of 
W and J~^W J . Then there exists a decomposition in space Vi of dimension 
2 with the Vi stable by rotations Rq. We can then write the variational 
equation under the following form 

t(-C2 + 2t + 2Et^)X + (-t + G'^)X = R-^^^AiReit)X i = l..n 

with y4j a 2 X 2 matrix (we can choose Ai diagonal after a basis change). 
Among the matrices Ai, there is one corresponding to the motion of the center 
of mass and this gives Ai = 0. There is also a matrix corresponding to the 
first integrals of the energy and angular momentum, and this corresponds to 
A2 = diag{2, —1). The other matrices do not have a priori special properties. 
Then the Galois group for the cases z = 1, 2 is C, and for the others, it is at 
most Sp{2). □ 

Theorem 14. Let V be the potential of the n body problem with positive 
masses mi, c an aligned Darboux point with multiplier —1 and W G M2„(C) 
the associated matrix (given by equation ( |T3l) ). The variational equation near 
a conic orbit is partially decoupled if and only if det(W) = 0. 

Proof. For an aligned Darboux point, we have the equalities f|T5|l . We pose 
V = {wi, W2)- If V is an eigenvector of W , then wi is an eigenvector of A and 
— \A with the same eigenvalue. Then det(y4) = 0. Conversely, if det(iy) = 0, 
then it exists an eigenvector w of eigenvalue of A, and then v = {w,w) is 
admissible. □ 
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5.2. The 3 body problem and some specific cases 

We already know that in all cases, the matrix W should have a double 
eigenvalue. Our approach will be the following. We search masses and Dar- 
boux points such that W has a double eigenvalue. Then for the corresponding 
eigenvector v, there are two possibilities. 

• Either Jv is also an eigenvector of W with the same eigenvalue. This 
corresponds to the case where the associated eigenspace is of dimension 
greater than 2. 

• Either v can be written v = {w,iw), and the matrix W is not diago- 
nalizable. 

For the aligned case, it is easier because we just have to look at the de- 
terminant. But in fact for the real ones, there are never zero eigenvalue if 



the Darboux point is real (this is due to the result of [12|), so we need to 
look at complex cases. But even there this constraint is much stronger than 
expected. We find the following theorem 

Theorem 15. Let V be the potential of the 3 body problem with positive 
masses mi,m2,m3 such that nii + m2 + = 1. Then V possess a Dar- 
boux point such that the variational equation near a conic orbit is partially 
decoupled if and only if 



{1711,1712,1713 




1 ^21 + V 126 + 42 V2T 1 
4 ^ 84 ' 2 ~ 

or permutation of these cases. 

Proof. Let us begin with aligned case. After renormalization, we can take 
c = (—1, 0, p) with p ^ 0,-1 and we have the Euler quintic equation 

L = {—mi — 7722) p^ + (— 3mi — 2m2) p'^ + (—3 mi — m2) 

(3m3 + m2) p^ + (3m3 + 2m2) p + m2 + ms = 

We search the eigenvalues of W , and we find that det{W) = if and only if 

2p2 + 3p + 2 = 
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After taking the resultant, we have 

Res{2p^ + 3p + 2, L, p) = 7m\ — 35mim2 — 35m2m3 + 56m^ + 63mim3 + bQm\ 

We want that this resultant vanish, and the only possibility for real positive 
masses is 

(mi,m2,m3) = 

We can permute the masses in the equation and this gives all the possibles 
permutations of this solution. But there is still a "complex order" and the 
corresponding potential is the following 

^ mim2 rriims 17121713 

qi -q2 qi- qs ?2 - ^3 

The Darboux point equation leads to 

L = (-mi - 1712) + (-3 mi - 2 m2) + (-3 mi + 2 ma - m2) 

(-2 mi + 3 ma + m2) + (3 ma + 2 m2) p + m2 + ma = 

The eigenvalues of W never vanish in this case. Let us look now at the 
Lagrange configuration. For complex coordinates, this corresponds to the 
case 

rf = 4 = rl 

where ri,r2,ra are the mutual distances between the bodies. We begin by 
the case ri = r2 = ra. We need a double eigenvalue and we find the condition 

3m2 — 3m2ma — 3mim2 + 3m3 — 3mima + Srnl — 

whose only solution is 

{mi, 1112,1713) = Q, ^, 

We check that the associated eigenspace of eigenvalue 1/2 is invariant by J, 

and it is the case. 

Let us look now at the complex cases. Among the 27—1 possibilities lots of 
them are in fact the same after dilatation permutation. After these reduc- 
tions, we find that there are only 3 essentially different cases 

(ri,r2,ra) = (1, l,j), (1, l,f), {l,j,f) j = e"^ 



(---) 
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The last one is also an aligned Darboux point (it is both Lagrange and 
Euler configuration), and so it already has been treated. First we search 
for masses such that W has a double eigenvalue. We find for (1,1, j) and 
(1, 1, j^) only one real positive solution, which is the last one of f lT6|) . This 
is the same for both Darboux points because they are conjugated. We look 
at the corresponding eigenspace (the double eigenvalue is 1/2), and we find 
that the eigenspace is only 1-dimensional. This is not enough for the case 
dim{V) > 2. In the case dim{V) = 1, we know that W should be non- 
diagonalizable. Moreover, the eigenvector should be written v v = {w,iw). 
We check these properties and they are satisfied. □ 

Remark 6. The last case of f|T6l) is very interesting for many reasons. We 
can study the variational equations and the structure of the equations is not 
so degenerated as in the other cases. Because of this, a more deeper analysis 



should be possible. For example, in [13|,[1J], another notion of partial inte 



grability is considered about the existence of a single additional first integral. 
For this last masses case, the two notions could probably be fused together 
to prove the non existence of a single additional first integral restricted to 
a single level of energy and angular momentum. This is because the varia- 
tional equation on the characteristic space associated to the eigenvalue 1/2 is 
simple enough to allow complete study, but is not trivial. Moreover, the fact 
that these masses do not possess any symmetry will avoid to consider special 
invariant sub manifold as the isosceles 3 body problem in, for example, the 
complete search of algebraic invariant manifold for the 3 body problem with 
these masses. 

Theorem 16. We consider V the potential of the n body problem in the 
plane with positive masses, and c a real Darboux point such that there exist 
a rotation 

Re, e i {kii, keZ} 

in the plane such that Rg sends the configuration on itself (conserving also the 
masses). Then there exists a double eigenvalue and the associated eigenspace 
is of dimension > 2. 

Proof. Let Rg be a rotation such that 6 ^ {kn, k E 1^} and that Re sends to 
configuration c on itself and conserve the masses. We pose W{c) as in flT3]) 
and we have then the identities 

W{Rec) = R_eW{c)Re W{Rec) = W{Pc) = p-^W{c)P 
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with P a permutation matrix (the rotation conserve the configuration and 
the masses of the bodies, but not the numeration of the bodies). Then 

W{c) = {ReP-'r'W{c)ReP-^ 

Let V be an eigenvector of W{c). Then RgP~^v is also an eigenvector with 
the same eigenvalue. We just have to prove it is not the same. We can write 
in a good basis 

f cos^/„ -sin^/„ \ Pa \ 

^ sin^/„ cos^/„ J \ P, J 

with P„ a permutation matrix. Wc have then that P and Re commute. We 
know that the rotation Rg is of finite order (because there are only a finite 
number of bodies and that the configuration is real). Then 6 = 27r/k with 
A; e N*, and k > 3. The matrix P is then also of order k. 

Let us consider the body number i with coordinates g^. We look at the 
orbit R^qi, j — 0...k — 1. This orbit contains either k elements or only one 
(and this case could only happen once, for a body placed on the center of 
mass) . We conclude that the permutation matrix should be of the following 
form 



(T 


... 


OA 








(T 









\ 





... 









Pa 




















... T 







or 











T 
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1 
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We conclude that the matrix RgP can be diagonalized in the form 

ReP ~ diag (e^^ e-^^ (e^(^■+^)^ . . . , e'^+'^', e'^^-'^', e'^^-'^').^^^_-^) 
or ReP ~ diag ((e^(^•+^)^ . . . , e^(^+^)^ e^(^-^)^ . . . , e^^^-'^') ^^^ ^_^) 

We suppose that the masses are positive and that the Darboux point is real, 
then all eigenvectors v of W{c) are real. Suppose that W{c) does not have 
any eigenspace of dimension > 2. Then all its eigenvectors are eigenvectors 
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of RqP. As RqP is real, if f is a real eigenvector of ReP-, then the associated 
eigenvalue is real and so the associated eigenvalue is ±1. This would mean 
that 

Sp{ReP) C {-1, 1} 
This is impossible because k >?>. □ 

5.3. The equal masses case 

Theorem 17. Let V he the potential of the n body problem in the plane with 
equal masses, c the Darhoux point given by the following 

f 2'K{l-l) \ . f 27,{l-l) \ . 

Ci = acos Ci+„ = asm t = 1 . . .n 

\ n J \ n J 

where a is such that the multiplier equals to —1. Let v be the vector given by 

/AnU - 1)\ fAnii - 1] , 

Vj = cos I f j+„ = sin I i = 1 . . .n 

n \ n 



Then (IT^ is satisfied with 



A = 2 



1 - ^os (^) sin (^) 



Proof. The proof is only a direct computation of the matrix W and then of 
Wv and the use of (lots of) trigonometric formulas. □ 

We can now eventually prove Theorem [2l 

Proof. Using Theorem [3, one just need to avoid specific values for A. We will 
then build a majoration and minoration for A given by formula f|T7|) . First 
of all, we remark that for n > 3 



'n-l N -1 



l-cos(^) \^sin(^) 



> 



Then A < 2. Let us prove now that A > 0. First we prove the following 
inequality 

sm{z) < g]0, 7r/2[U]7r/2, 7i[ 

Sill \^) 

27 



and we compute the formula 



/7r\ "—1 / 



/ N sin f — 

l-cos(^) ^ ^ " 



Using both of them, this gives for n > 3 

{-) ^ 1 
1 - (^) ^ Sin (=f ) 

So we get that A > 0. Using the integrabihty table of Theorem (TJ there are 
no exceptional values in ]0, 2[. □ 

The case C^H = is special. In this case, we have either C = or H = 
(or both). The case H = corresponds to parabohc orbits. These orbits are 
used by Tsygvintsev in Q, and he solves the case for 3 bodies with equal 
masses (he studies also existence of a single additional first integral that we 
do not consider). In the case of C = 0, the problem is solved by Morales, 
Simon in [lij for the n equal masses. Our reasoning is also valid for all these 
cases. 

Remark 7. The only left case is H = C = 0. Here, the variational equation 
is always integrable and in fact it is always the case at all orders. This is linked 
to the fact that we can reduce the system using homogeneity and rotation, 
allowing to diminish the dimension of 4. We obtain then a "direction" field 
(we lose notion of time after reduction, but not the integrabihty notions) on a 
manifold of dimension 4n — 8. This, however, destroy Hamiltonian structure, 
and moreover, the Darboux points correspond now to fixed points of this 
field. One would need a new particular orbit (explicit) to apply Morales 
Ramis method, but no such orbit is known. 

I wish to thank Andrzej Maciejewski and Maria Przybylska who, maybe 
without even noticing it, help me to enhance this work by their remarks, 
discussions and invitations in Zielona Gora. 
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